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In this paper we study static spherically symmetric Einstein-Vlasov shells, made 
up of equal mass particles, where the angular momentum L of particles takes values 
only on a discrete finite set. We consider first the case where there is only one value 
of L, and prove their existence by constructing explicit examples. Shells with either 
■ hollow or black hole interiors have finite thickness. Of particular interest is the thin 

Q_|. shell limit of these systems and we study its properties using both numerical and 

analytic arguments to compare with known results. The general case of a set of 
values of L is also considered and the particular case where L takes only two values 
is analyzed, and compared with the corresponding thin shell limit already given in 
the literature, finding good agreement in all cases. 
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I. INTRODUCTION 



^ \ Although sometimes sidestepped, it is a general requirement in studying non vacuum 

spacetimes in general relativity that the energy momentum tensor, that is the matter (field) 
contents, should have a clear, although possibly highly idealized, physical interpretation. 
ON ! Among these choices the case where matter is described as a large ensemble of particles that 
interact only through the gravitational field that they themselves, at least partially, create, 
is of particular interest, both because of their usefulness in modeling physical systems such 
as star o galaxy clusters, and of the possibility of a relatively detailed analysis, at least 
in some restricted cases. As usual in theoretical treatments, one starts imposing as many 
restrictions as compatible with the central idea, and then tries to generalize from these 
cases. In this respect, the restriction to static spherically symmetric systems provides an 
important simplification, although even with this restriction the problem is far from trivial, 
and further restrictions have been imposed in order to make significant advances. One of the 
first concrete examples is that provided by the Einstein model (H, where the particles are 
restricted to move on circular orbits. This model is static, and it is not easy to generalize 
as such to include a dynamical evolution of the system. This generalization can, however, 
be achieved if the particle world lines are restricted to a shell of vanishing thickness( "thin 
shell"), as considered by Evans in (2|. The analysis in 0], although, motivated by the 
Einstein model, considers only shells where all the component particles have the same value 
of their (conserved) angular momentum. In a recent study [l| of the dynamics of spherically 
symmetric thin shells of counter rotating particles, of which j2| is an example, it was found 
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that the analysis can be extended to shells where the particles have angular momenta that 
take values on a discrete (but possibly also continuous) set, and is not restricted to a single 
value. It was also found that in the non trivial thin shell limit of a thick Einstein shell the 
angular momentum of the particles acquires a unique continuous distribution, and, therefore, 
the models in [2J and [lj are not approximations to an Einstein model. A relevant question 
then is what, if any, are the (thick) shells that are approximated by those in ^ and In this 
paper we look for an answer to this question by considering a generalization of the Einstein 
model where instead of circular orbits we impose, at first, the restriction to a single value 
of the angular momentum. The particle contents is described by a distribution function 
/ in phase space, and, because of the assumption of interaction only through the mean 
gravitational field, / must satisfy the Einstein- Vlasov equations 0. In the next Section we 
set up the problem and show that it leads to a well defined set of equations. In Section III 
we set up and analyze a particular model, obtaining expansions for the metric functions at 
the boundary of the support of /, appropriate for numerical analysis. Further properties are 
analyzed in Section IV, where we show that all these shells have finite thickness. Section 
V contains numerical results for a generic example. The "thin shell" limit is considered in 
Section VI, both through analytic arguments and a concrete numerical example, with the 
results showing total agreement with the thin shell results of [lj. A further comparison with 
[l| is carried out in Section VII, where the stability of a shell approaching the thin shell limit 
is considered. The generalization to more than one value of L is given in Section VIII, where 
we find that particles with different values of L may be distributed on shells that overlap 
completely, or do so only partially or not at all. Numerical examples and a comparisons 
with [lj are finally developed in Section IX. Some comment and conclusions are given in 
Section X. 



II. THE STATIC SPHERICALLY SYMMETRIC EINSTEIN- VLASOV SYSTEM 

The metric for a static spherically symmetric spacetime may be written in the form, 

ds 2 = -B(r)dt 2 + A(r)dr 2 + h(r) 2 dQ 2 (1) 

where dfl 2 = d6 2 + sin 2 9d(f> 2 is the line element on the unit sphere and r > 0. 

For a static, spherically symmetric system, the matter contents, in this case equal mass 
collisionless particles, is microscopically described by a distribution function F(r, p- 7 ), where 
p> = (p r ,p e ,p^) are the components of the particle momentum, taken per unit mass. Then, 
as a consequence of the assumption that the particles move along geodesies of the space 
time metric, the distribution function F satisfies the Vlasov equation, which, in this case, 
takes the form, 

p r d r F-Ti b p a p b d p] F = (2) 

where a,b correspond to (t,r,9,<p). It is understood in (j2j) that p l is to be computed using 
9abP a p b = —1, so as to satisfy the "mass shell restriction" fj, — 1, where \i is the particles 
mass. Therefore, in what follows we set, 



JB(r 



= \/l + A{r){p r ) 2 + h(r) 2 \{p e ) 2 + sin 2 6{p^) 2 ] (3) 



We also notice that p t = dt/dr, where r is proper time along the particle's world line. 
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The Einstein equations for the system are, 



Gab '■= R a b — -^Rg a b — %nT ab (4) 
with the energy momentum tensor given by, 

T a b = ~ / Fp a p b \g\ i (5) 

J Vt 

where g is the determinant of g a b, and p a = g a bP h ■ Equations (T5]), (J4j) and ([5]) define the 
Einstein- Vlasov system restricted to a static spherically symmetric space time, with the 
metric written in the form ([I]). 

The assumption that the metric is static and spherically symmetric implies conservation 
of the particle's energy, 

E = B(ry 

= v 7 ^) y 1 + A(r)(p r ) 2 + h(r) 2 [{p 9 ) 2 + sin 2 0{p*) 2 ] (6) 
and of the square of its angular momentum per unit mass, 

L 2 = h(r) 4 Up 6 ) 2 + sin 2 8 (p^) 2 ] (7) 
It is easy to check that the Ansatz, 

F(r,p>) = ®(E,L 2 ) (8) 

where E, and L 2 are the functions of r, and pP given by (JHE])? solves the Vlasov equation 
for an arbitrary function $. 

To construct and solve explicit models based on ([8]) for the metric ([1]), it is convenient to 
change integration variables in ([5]). We set, 



1 

h{r) 



P = T7Z^ Lcos X 



V* = h( VT - n Ls ™X (9) 
h{ry smS 



and write (El) in the form, 



/^( r ) fchtv r2\ LdLdxdp r 

T - (r) = Wfiw)l ( ' )PaPb — ? — (10) 



where we should set, 



V 1 = \li{l + (P r ) 2 A+^) (11) 



B V h 2 J 
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Andreasson and Rein have explored the properties of models where $ takes the form, 

®(E,L 2 ) = ( f ) (E/E )(L 2 -Liy (12) 

In this note we consider a different type of models, based on the Ansatz, 

L 2 ) = F(E) Q(E - E) 5{L - L ) (13) 

where Q(x) is the Heaviside step function, and 8(x) is Dirac's 5; namely, we assume that L 
takes only the single value L , and that there is an upper bound on E, given by E . F(E) 
is assumed to be a smooth function of E. We then have, 

AttT r\\/~A fPmax I 

T t l = - fe ° 3 V J F{E)^Jh? + L 2 + {prfh?Adf 

Tr = 47rWg ri™ F(E)(f) 2 d r 
h J ^h 2 + L 2 + (p^h 2 A P 
e = 2n(L ) 3 VA r*~ F(E) 



h 3 Jo ^h 2 + L 2 + (p r ) 2 h 2 A 
T/ = T e e (14) 



where p r max depends on r and is given by, 



fmax ~ \j A(r)\] B(r) 1 h(rf (15) 

if Eq > B(r)(l + Ll/h 2 ), and p^ ax = otherwise. This simply states the fact that T ab ^ 
only in those regions where a (test) particle with energy E and angular momentum L Q can 
actually move. 



III. PARTICULAR MODELS 

We may now use the previous results to construct simple models and analyze their inter- 
pretation for a range of possible parameters. This analysis may be carried out in a number 
of ways. Here we choose the following; we first use a gauge freedom in the metric ([TJ) to set, 

h(r) = r 

A(r) = 1/(1 - 2m(r)/r) (16) 



Then, from the Einstein equations and the form (JMJ) of T a b, we find two independent equa- 
tions for m(r) and B(r), 

dm 2 
— — = 47rr p[r) 
dr 

dB 2B(r)(m(r) + 4:7rr 3 p(r)) 
dr r(r — 2m(r)) 

where, p(r) = — T t * is the energy density, and p(r) = T r r is the radial pressure, given by (|14p . 
with h(r) = r. There is also an equation for pr(r) = Tg e , but, as can be checked, this is not 
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independent of (JTTj) . Equations (fT7|) are deceivingly simple, because the explicit dependence 
of p and p on m and £> is in general quite complicated. Here we consider a simple example 
and propose a method for constructing the solutions, that is illustrated by the example. It 
can be seen that some simplification is attained if we choose, 




F(E) = QiE = QUlB[l + (p r ) 2 A + ^ ) 



where Q\ > is a constant. With this choice we may perform the integrals in ( Fl4|) explicitly 
and, after some simplifications, we get, 

dm Q 2 [2(Lg + r 2 )B + r 2 E 2 } ygg| - (Lg + gg 

rffi = 2mB + 2Q 2 [r 2 L 2 - (Lg + r 2 )Lf /2 
dr r(r — 2m) r 3 (r — 2m) 

where Q 2 = 167r 2 L Qi/3. We also find, 



p(r) 
p(r) 



Q 2 [2(Lg + r 2 )B + r 2 E 2 ] y/r 2 L 2 - (Lg + r 2 )L? 
4irr 5 B 

Q 2 [r 2 E 2 -(L 2 + r 2 )Bf 2 
4nBr 5 



p T {r) = ^° VI a °_~ ^ (21) 



3Q 2 L 2 v /r*EZ-(L'i ) +-^jB 
87rr 5 

Considering ( [1911201) . we find that it is a simple, but rather difficult to handle, system of 
equations for m{r) and B(r). We have not found closed (analytical) solutions for the system, 
and, therefore, we must resort to numerical methods. The application of these methods 
requires, however, considering and solving several subtleties inherent in the system. As 
indicated above, in all these equations ( f)19ll20p and ([2~T]) ). the terms involving Q 2 should be 
set equal to zero if r 2 E^ < (Lg + r 2 )B. We notice that for Qi = Q 2 = we have m(r) = M 
with M = constant, and B[r) = B (l — 2M/r), where B is also a constant, corresponding 
to the standard Schwarzschild solution. For a shell type solution, these solutions correspond 
to the inner and outer regions, to be matched to the region where T a b 7^ 0. When Qi ^ 0, 
since we must have B(r) > 0, we must also have dm/dr > 0, but, even though p > we 
might end up with m(r) < 0, and still have all equations satisfied. 

For shell like solutions, either with an empty interior or with a central mass (black hole), 
a further difficulty can be seen considering that there should exist an "allowed region" where 
t 2 Eq > (Lg + r 2 )B, with r taking values in the interval T{ < r < r Q , where and r Q are, 
respectively, the inner and outer radii of the shell. We must impose continuity in both B(r) 
and m(r) to avoid S functions in T a t,. This implies that t 2 Eq — (Lg + r 2 )B is continuous 
in Ti < r < r Q and approaches continuously the value zero at the boundaries. Therefore, 
both dm/dr and dB/dr are also continuous inside and at the boundaries of this interval, 
and actually we have dm/dr\ r = r . = 0. We also find that d 2 B/dr 2 should be continuous, but 
d 2 m/dr 2 must be singular, and this makes the construction of numerical solutions where we 
try to fix from the beginning the values of T{ and r D rather difficult. Nevertheless, the above 
analysis indicates that for r > r^, but r ~ r iy we should have, 

B{r) = B + B 1 (r- n) + B 2 {r - r,) 2 + TZ B 

m{r) = M l +U m (22) 
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where B Q , Bi, and Mi are constants and TZb and TZ m are functions of r that vanish respec- 
tively faster than (r — Tj) 2 and (r — Tj) for r — r 4 — > + . It is straightforward to extend this 
analysis to higher order by imposing consistency between the right and left hand sides of 
ffigpOl) as r - n -> 0+. We find, 

B(r) = B + Bi(r- n) + 5 2 (r - r,) 2 + £ 3 (r - r;) 5 /2 + TZ B 

m(r) = Mi + M 2 (r - ri ) 3/2 + M 3 (r - r,i) 5/2 + TZ rn (23) 

where B 2 , B 3 , M 2 and M 3 are constants, and 7^.^ an d lZ m stand for higher order terms. The 
constants appearing in (123]) are not independent. They may be written, e.g., in terms of r*, 
Mi, Lq, Eq and Q 2 . We notice that Mi is the Schwarzschild mass for the region inside the 
shell (r < rj). Moreover, the system (I19II20I) is invariant under the the rescaling B —>■ \B, 
E —> y/XE , and Q 2 — > Q 2 / V^X. 

The condition that r = corresponds to the inner boundary of the shell implies, 

Bo = ^ (24) 



Similarly, we find, 

Bi 
Bo 



2MiEln 



(r t -2Mi)(rf + L 2 



0) 



2 



MiEq 



(r i -2M 1 )(r 2 + L 2 ) 



M 2 = Q2 E M±M^ zffizM (25) 

The explicit expressions for B3 and M3 are also easily obtained but are rather long and 
will not be included here, although they were used in the numerical computations described 
below. 



IV. SOME GENERAL PROPERTIES 

An interesting question regarding the model of the previous Section is related to the 
possible values that the thickness of the shells can attain. This may be analyzed by con- 
sidering the limit of solutions of the system (I19|l20p as r — > 00, under the restrictions that 
p 7^ 0, and r > 2m(r). The first, according to (I19ll20p and (l2Tj) . implies B(r) < Eq, and 
Bq = lim^oo B (r) < Eq. We remark that dB/dr > 0. Therefore, B(r) must approach Bo 
monotonically from below. Consider first the case Bq < Eq. Replacing in the first equation 
in (I19|I20I) . for large r we find that dm/dr approaches a constant value, and, therefore, m(r) 
grows linearly with r. But then, replacing in the second equation in (I19|I20I) . we find that 
dB/dr decreases as 1/r, leading to a logarithmic growth in B(r), incompatible with the 
assumed conditions. Therefore, any possible solution should have Bq = Eq. Then, for large 
r we should have, 

B(r) = -^L - Bi(r) (26) 
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with B\(r) — > 0, as r — > oo. Replacing now in fll9|20p . to leading order we find 

dm 



- 3Q 2 V£i (27) 

and this implies that m(r)/r — > as r — > oo. Then, using again fll9|20jl . we should have, 

dBi 2£ 3 Q 2 



(28) 



and this implies B\(r) ~ —2EqQ 2 ln(r), which contradicts the assumption B\(r) — >• 0. Thus 
we conclude that the equation, 

£ V - (Lq + r 2 )5(r) = (29) 

must always be satisfied for finite r, and, therefore, all shells constructed in accordance with 
the prescription [Hfy have finite mass and finite thickness. 

We nevertheless believe that this results is more general, and applies to all shells satisfying 
the Ansatz (1131) . although we do not have a complete proof of this statement. 



V. NUMERICAL RESULTS 

As indicated, we do not have closed solutions of the equations for B(r) and m(r), even 
for the simple model of the previous Section. Nevertheless, since (I19) 20p is a first order ODE 
system, we can apply numerical methods to analyze it. We may use the expansions (1231) . 
(disregarding the terms in TZb and TZ m ), to obtain appropriate initial values for B{r) and 
m(r), for r close to r i; in the non trivial region r > r^. 

We may illustrate this point with a particular example. We take rj = 7.0, Mi = 1.0, 
L = 4, E = 1 and Q 2 = 0.1. Using these values and (I23p (truncated as indicated above), 
we find 5(7.0001) = 0.7538..., and m(7.0001) = 1.0000... (actually, the computations were 
carried out to 30 digits, using a Runge-Kutta integrator). The numerical results are plotted 
in Figure 1 and Figure 2. 



VI. THIN SHELLS 

One of the motivations for studying the type of shells considered in this paper is the 
possible existence of a non trivial "thin shell" limit, where the thickness of the shell goes to 
zero, with the restriction to a single or a finite set of values of L, and how this limit compares 
with the thin shells considered in [2| and [l|. We remark that the existence of thin shell 
limits of Einstein- Vlasov systems has already been analyzed in the literature Q. Here we 
are interested not only in the existence of this limit for our particular models, but especially 
in the limiting values of the parameters characterizing our shells. Since this type of analysis 
is not immediately included in, e.g, Q, we consider it relevant to provide an explicit proof 
of the properties of our models in the thin shell limit. 

We first recall that for a static thin shell constructed according to Evans' prescriptions 
[2|, we have the following relation between the radius R, inner (Mi) and outer (M 2 ) mass, 
and angular momentum L of the particles, 

~ 2 = R- ^/R-2M^R-2M 2 
3VR-2M IV /R-2M 2 - R 1 ' 
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FIG. 1: Plots of B(r) (solid), m(r) (dotted), and 30(dm/dr) (dashed), as functions of r for rj = 7.0, 
Mi = 1.0, L = 4, E = 1 and Q 2 = 0.1. The end point of the plot is at r Q = 14.1098..., with 
B(r a ) = 0.9256... and M 2 = m(r ) = 1.2306... 
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FIG. 2: Plots of p{r) (solid, higher curve), 4p^(r) (dashed, middle curve), and 40p(r) (dotted, 
lower curve), as functions of r for rj = 7.0, M\ = 1.0, Lq = 4, Eq = 1 and Qi = 0.1. The end point 
of the plot is at r Q = 14.1098..., with B(r a ) = 0.9256... and M 2 = m(r ) = 1.2306... 
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We may now prove that the non trivial thin shell limits of the shells constructed according 
to the prescription (Tl8|) effectively coincide with the Evans shells of reference [2| as follows. 
We first take the r derivative of (120]) . and then use (Tl9|) and (1201 to obtain, 



d 2 B 
dr 2 



r{2E 2 r 2 



rl B 

L 2 )B)^f- + 2(AE 2 r 2 + (2L 2 - r 2 )B) 
dr 



dm 



r(E 2 r 2 + 2(L 2 + r 2 )B)(r - 2m) 



dr 



r 2 (r — 2m) 



9 \ dB 

r(r — 2m)—— 

dr 



mB 



(31) 



m r,- 



Next let rj and r be, respectively, the inner and outer radii of the shell, and Mi 
and M 2 = m(r ) the corresponding masses inside and outside the shell. Then, for r, < r < r c 
we have dB/dr > 0, and, 



dB 
dr 



2M 1 B(r i 



ri nin - 2M l ) 

The idea now is to use the fact that, 



dB 
dr 



2M 2 B(r Q ) 
r (r o - 2M 2 ) 



r ° d 2 B , dB 
dr = —— 



dr 2 

r ° dm , 
——dr 
dr 



dr 
M 2 



dB_ 

dr 



2M 2 B{r ) 



2M 1 B(r l 



r (r - 2M 2 ) nfc - 2M l ) 



Mi 



(32) 



(33) 



On this account we rewrite (l3T|) in the form, 

1 dm r(E 2 r 2 + 2B(L 2 + r 2 )) 



d 2 B 



r — 2m dr 



r(2E 2 



L 2 )B) 



dB 

dr 



2(AE 2 r 2 + (2L 2 - r 2 )B) 



dr^ 



A(E 2 r 2 + 2B(L 2 + r 2 )) [r(r -2m)^§- mB] 



r(r — 2m) 



r(2E 2 + (r 2 + L 2 )B) 



dB 
dr 



2(AE 2 r 2 + (2L 2 - r 2 )B) 



(34) 



and integrate both sides from r = Tj to r 



But now we notice that while both m(r) 



and dB/dr are rapidly changing but bounded in < r < r G , the change of S(r), and r in 
that interval is only of order r D — r^. We may then choose a point r = R, with < R < r a , 
and set B(r) = B(R) = B , and r = R, except in the arguments of m(r), dm/dr, dB(r)/dr 
and d 2 B(r)/dr 2 , in (I3"l"|) . as this introduces errors at most of order r Q — r i; in the factors of 
dm/dr, and d 2 B(r)/dr 2 , and in the last term in the right hand side of (1341) . Similarly, we 
may set E 2 = [L 2 + R 2 )B Q /R 2 in (I3"41) . to obtain, up to terms of order r Q — r i; 



dm 



R(L 2 + R 2 )) 



d 2 B 



R — 2m{r) dr 



dB rlr 2 

R(R 2 + L 2 ))—— + 2(i? 2 + 2L 2 )B ar 
dr 

4((L 2 + R 2 )) [R(R - 2m) f - mB ] 



+ 



R(R - 2m) 



d B 

R(R 2 + L 2 )— + 2(R 2 + 2L 2 )B 



(35) 



10 




~7~ 7.002 7.004 7.006 7.008 7.01 7.012 7.014 



FIG. 3: Plots of m(r) (dotted, higher curve), B(r) (solid, middle curve), and Wp (dashed, lower 
curve), as functions of r for = 7.0, Mi = 1.0, Lq = 4, Eq = 1 and Q 2 = 800. The end point of 
the plot is at r Q = 7.0155..., with B(r Q ) = 0.7546... and M 2 = m(r ) = 1.2999... 



and, again, we notice that the last term on the right of (|35|) gives a contribution of order 
r G — rj. We then conclude that, 

lim f — rr~^dr = lim [ —-75 — — — r^r^ r (36) 

R - 2m(r) dr J n + d_B + + ^ dr* 

dr 



The integration of the terms in dm/dr, and d 2 B(r)/dr 2 , is now straightforward. We use 
next ( |32|) and the fact that in this limit B(ri) = Bq = B(r ) and rj = R = r Q , to obtain, 

(37) 



y/R - 2Mi _{R- 2M 1 ){R 3 + 2RL 2 - 3M 2 L 2 - M 2 R 2 ) 
y/R - 2M 2 ~ (R - 2M 2 )(R 3 + 2RL 2 - "iM x L 2 - M X R 2 ) 

Solving this equation for L 2 , we finally find, 



-2 



R 2 (R - y/R-2M 1 yJR-2M 2 ) 
3VR ~ 2M[y/R -2M 2 -R 



(38) 



which is, precisely, the relation satisfied by the parameters of the shells of ( 1301 ). 

We can also check this result, and, in turn, the accuracy of numerical codes, by directly 
considering initial data for the numerical integration that effectively lead to shells where the 
thickness is a small fraction of the radius. 

A particular example is given in Figure 3, where the values of the initial data is also 
indicated. We can see that the mass increases by about 30 percent, while the thickness of 
the shell is less than 0.3 percent of the shell radius. We can check that these results are in 
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agreement with (1301) . Solving this for M 2 , 

R \2(R - 3M 1 )R 2 L 2 + (4R - VM^Ll - M X R A 
Mo = —± 



(R-2M 1 )(R 2 + 3L 2 ) 2 



(39) 



and replacing L = 4, Mi = 1, and R = 7, we find M 2 = 1.2997... in very good agreement 
with the numerical results quoted in Figure 3. 

VII. STABILITY ANALYSIS 

Consider a shell approaching the thin shell limit. We restrict to the case of vanishing 
inner mass Mi = 0. If r« and r Q are, respectively the inner and outer boundaries of the shell, 
the matching conditions in the absence of singular shells at and r Q imply that both m(r) 
and B(r) and their first derivatives should be continuous at both r = r« and r = r Q . The 
form (JTJ) of the metric with the choice (fT6l) imply that for < r < r« we have, 

ds 2 = -Bidt 2 + dr 2 + r 2 dU 2 (40) 

where B^ > is a constant. Then, at r = Tj we have, 

r 2 E 2 

Bi = i£S (41) 
For r > r„, the metric takes the form, 



,.o (. 2M^ 



r 



ds * = - Bo [ t _L_ ) dt 2 + ( 1 - ^P- ) dr 2 + rW (42) 



where B > is another constant satisfying, 



r 2 E 2 



L 2 + r 2 o 

Consider now a test particle moving along a geodesic of the shell space time, with 4- 
velocity C/ M (r) = (dt/dr, dr/dr, dtO/dr, dtfi/dr), with U^U^ = —1. Without loss of generality 
we may choose dO/dr = 0. Then we have the constants of the motion: 

E = B{r)— , L = r 2 ^ (44) 
dr dr 

and the normalization of implies, 

Therefore, for r > r Q we have, 

dr\ 2 E 2 ( 2Mo\ ( L 2 \ ( 2M- 
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and the particle radial acceleration is given by, 

d 2 r L 2 (3L 2 + r 2 )M 2 



dr 2 r 3 r 4 



(47) 



If we assume now that the shell is close to the thin shell limit, with angular momentum L$, 
radius R, (with r, < R < r Q ), and mass M 2 , then we should have, 



^H- VR-2M 2 )R 2 



Then, for L ~ L , and r ~ R ~ r Q we find, 

d 2 r - 2M 2 ( v / i? - v^R - 2M 2 



Lq = vv ; v (48) 

3Vi? -2M 2 -VR 



< (49) 



rfr 2 i?3/2 (3^ _ 2M 2 - v^R) 

Similarly, in the region < r < r^, for r ~ ~ R we find 

d 2 r 2(^fR- y/R-2M 2 



dr 2 R(3y/R - 2M 2 - v 7 ^) 



> (50) 



and, therefore, all these shells are stable under "single particle evaporation", in total agree- 
ment with the results of [ill. 



A related problem is that of the dynamical stability of the shell as a whole, as was also 
analyzed in |l|. There, the shells considered where "thin", and therefore, the motion was 
described by ordinary differential equations for the shell radius R as a function, e.g., of 
proper time on the shell, which allowed for a significant simplification of the treatment of 
the small departures from the equilibrium configurations. Unfortunately, the corresponding 
equations of motion for the shells considered here would be considerably more complicated 
and their analysis completely outside the scope of the present research. We, nevertheless, 
expect that such treatment, if appropriately carried out, would also agree with the results 
found in [l[ as the "thin shell" limit is approached. 



VIII. SHELLS WITH TWO OR MORE VALUES OF THE ANGULAR 

MOMENTUM 

It is rather simple to extend the analysis of the previous Sections to the case where the 
angular momentum of the particles takes on a discrete, finite set of values. Instead of fTTBl . 
we have, 

L 2 ) = F t {E) Q(Ei - E) 5{L - U) (51) 

i 

where the functions > are arbitrary, with i = 1,2, ...,N, and iV finite. We will restrict 
to the case of two separate values, (N=2), since, as will be clear from the treatment, the 
extension to a larger number of components is straightforward. We are actually interested 
in the behaviour of these shells as they approach a common thin shell limit. Therefore, we 
will further simplify our Ansatz to the form, 




Fi(E) = QiE = Q iX B 1 + (p^ A + -i (52) 
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where Qi > are constants. With this choice we may perform the integrals in (fT4"I) explicitly 
and, after some simplifications, we get, 

dm d [2{L\ + r 2 )B + r 2 E 2 } ^r 2 E\ - {L\ + r^B 

dr r 3 B 

C 2 [2(Lj + r 2 )B + r 2 E 2 } ^r 2 E 2 -(L 2 + r 2 )B 
+ ^Fg 

dB 2mB 



dr r(r — 2m) 

2d \r 2 E\ - [L\ + r 2 )Bf /2 2C 2 [r 2 E 2 - (L 2 + r 2 )B} 3/2 
r 3 (r — 2m) r 3 (r — 2m) 

where Ci = 16ii 2 LiQi/3, % = 1,2. We also find, 



(53) 



p{r) 



Ci [2(L 2 + r 2 )B + r 2 E\\ ^r 2 E 2 - (L\ + r 2 )5 
47rr 5 S 



p(r) 



C 2 [2(L 2 + r 2 )5 + r 2 E 2 } ^r 2 E 2 - [L\ + r 2 )B 
4irr 5 B 

Ci [r 2 E 2 - (L 2 + r 2 )£?] 3/2 C 2 [r 2 E 2 - (L 2 + r 2 )ff] 3/2 
4 n Br 5 + 47r.Br 5 



3C 1 L 2 v /r 2 E 1 2 -(L 2 + r 2 )i? 3C 2 L 2 ^r 2 E 2 - (L 2 + r 2 )g 
87rr° 87rr° 

It is clear that we recover the results of the previous Sections if we set either C\ or C 2 
equal to zero. 

It will be convenient to define separate contributions to the density, p\ and p 2 , for the 
particles with L\ and L 2 . 



Ci [2(L 2 + r 2 )B + r 2 E 2 } ^r 2 E\ - (L\ + r 2 )B 
Pl{T} ~ Anr 5 B 



C 2 [2(L 2 + r 2 )B + r 2 E 2 } ^r 2 E 2 - [L\ + r 2 )B 
P2[r) ~ 4^B (55) 

Then, provided the integrations cover the supports of both pi and p 2 , we have, 

M 2 - Mi = AM = Ami + Am 2 (56) 

where 

Ami = J 4:iir 2 pidr , Am 2 = J ^ixr 2 p 2 dr (57) 

These expressions may be considered as the contributions to the mass from each class of 
particles. This will be used in the next section to compare numerical results with the thin 
shell limit. 

Equations fl53l) may be numerically solved for appropriate values of the constants Cj, E i: 
and and initial values, i.e., for some r, of m(r) and B(r). We remark that, as in the 
previous Sections, it is understood in fl53|) . (and also in (1541) ). that both ^Jr 2 E 2 — [L 2 + r 2 )B 
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and [r 2 E 2 — [L 2 + r 2 )B] 3 ^ 2 must be set equal to zero for r 2 E 2 < (L 2 + r 2 )B. In this general 
case, it is clear that the shells (where by a "shell" we mean here the set of particles having 
the same angular momentum Li) may be completely separated or they may overlap only 
partially. We are particularly interested in the limit of a common thin shell for the chosen 
values of Li. One way of ensuring that at least one of the shells completely overlaps the 
other is the following. We choose an inner mass Mi, and an inner radius r^. This implies 
m(ri) = M\, while B(ri) = Bi is arbitrary. If we choose now arbitrary values for Li and L 2 , 
the density p will vanish at r = if we choose, 



El 



B,{L\ + r 2 ) 

El = B 'K +r ?> (58) 



Actually we also need to impose, 



Ll > 



r 2 M x 
n - 3Mi 



T 2 Mi 

12 > W^M, (59) 

to make sure that r = ri is the inner and not the outer boundary of the shells. We shall 
assume from now on that L2 > L\. Since, using the same arguments as in the single shell 
case, the shells have finite extension, it follows that one of the shells will be completely 
contained in the other. 

In the next Section we display some numerical results, both for thick shells that overlap 
partially, and for shells approaching the thin shell limit. We again find that the limit is 
associated to large values of the C iy and that the parameters describing the shells approach 
the thin shell values found in hi. 



IX. NUMERICAL RESULTS FOR TWO COMPONENT SYSTEMS 

As a first example we take Li = 5, L 2 = 6.5, = 7, and Mi = 1. We also set B{rj) = 1, 
for simplicity. Then, from (1581) . we set E\ = 1.510..., E 2 = 1.862... Finally, we choose 
Ci = 5, C2 = 3, and carry out the numerical integration. The results obtained indicate that 
the particles with Li — 5 are contained in the region 7.0 < r < 7.2682.., while for L2 = 6.5 
the corresponding range is 7.0 < r < 7.4764.... The resulting value of the external mass in 
M 2 = 2.1338.... In Figure 4 we display the total density p as a function of r (solid curve), 
as well as the contributions pi and p 2 to the density from the particles with respectively 
Li = 5 (dashed curve) and with L 2 = 6.5 (dotted curve). 

As an illustration of the approach to the thin shell configuration we considered again 
the previous values Li = 5, L 2 = 6.5, r { = 7, Mi = 1, -B(rj) = 1, but choose Ci = 800, 
C 2 = 240, and carried out the numerical integration. Figure 5 displays the functions p(r) 
(solid curve), pi(r) (dashed curve), and p 2 (r) (dotted curve). We see that now the shell 
extends only to the region rj = 7.0 < r < r Q = 7.020..., i.e., its thickness is less than one 
percent of its radius. The mass, on the other hand, increases by roughly a factor of two, 
since M 2 = m(r ) = 2.022.... 
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FIG. 4: Plots of p(r) (solid), pi(r) (dotted), and p 2 ( r ) (dashed), as functions of r for = 7.0, 
Mi = 1.0, Li = 5, L 2 = 6.5, ^ = 1.51..., E% = 1.86.. and C x = 5, C 2 = 3. The end point of the 
plot is at r Q = 7.4764..., with B(r ) = 1.0606... and M 2 = m{r ) = 2.1338... 




FIG. 5: Plots of p(r) (solid), p\{r) (dotted), and p2(r) (dashed), as functions of r for = 7.0, 
Mi = 1.0, Li = 5, L 2 = 6.5, E\ = 1.51.., E\ = 1.86.. and d = 800, C 2 = 240. The end point of 
the plot is at r Q = 7.020..., with B(r ) = 1.0026... and M 2 = m(r ) = 2.022... 
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We may compare these results with those of the thin shell limit of |l| as follows. It can be 
seen from (29) and (31) in [l| that for a thin shell of radius R, inner mass M\ and outer mass 
M 2 , with two components with angular momenta L\, and L 2 , the ratio of the contributions 
of each component to the total mass is given by, 

Am 2 = (L*-L*)(R* + LZ) 
Am, (L2-Zg)(iP + L?) 

where L is given by (I3"U|) . The numerical integration gives Ami = 0.4545..., Am 2 = 
0.5708.... If now take R = 7, and solve f[6"0"j) for L we find, L = 5.8079..., while replacement 
in (130]) gives Lq = 5.8386..., which we consider as a good agreement, with a discrepancy of 
the order of the ratio of thickness to radius. 



X. FINAL COMMENTS AND CONCLUSIONS 

The general conclusion from this work is that one can effectively construct a wide variety 
of models satisfying the restriction that L takes only a finite set of values, and that they 
do seem to contain the models used in [l| as appropriate thin shell limits. We remark also 
that the starting point for our construction is a variant of the Ansatz used in Q, where 
/ is factored in an E (the particle energy) and an L dependent terms. The possibility of 
multi-peaked structure in the case of more than one value of L obtained here is also in 
correspondence with the general results obtained in 0]. 
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